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We present an approach to calculate the atomic bulk solubility in binary alloys based on the 
statistical-thermodynamic theory of dilute lattice gas. The model considers all the appropriate 
ground states of the alloy and results in a simple Arrhenius-type temperature dependence determined 
by a "low-solubility formation enthalpy". This quantity, directly obtainable from first-principle 
calculations, is defined as the composition derivative of the compound formation enthalpy with 
respect to nearby ground states. We apply the framework and calculate the solubility of the A specie 
in A-Ti alloys (A=Ag,Au,Cd,Co,Cr,Ir,W,Zn). In addition to determining unknown low-temperature 
ground states for the eight alloys, we find qualitative agreements with solubility experimental results. 
The presented formalism, correct in the low- solubility limit, should be considered as an appropriate 
starting point for determining if more computationally expensive formalisms are otherwise needed. 

PACS numbers: 



I. INTRODUCTION 

High-throughput ah initio methods, capable of predict- 
ing properties of an ample set of materials from quantum 
mechanics calculations, are becoming important tools 
for scientists working in rational materials development. 
These methods allow researchers to correlate between dif- 
ferent systems and to observe trends convergmgtoward 
the predictions of new materials [J, H, d, 0,11, HI- The 
main difference between the "several-calculations" and 
the "high-throughput" philosophies is that the latter re- 
quires rapid estimations of materials properties so that 
the correlations between systems, even if roughly char- 
acterized, become the target information instead of the 
throughout and accurate description of a small subset. 
Clearly, to analyze the extensive amount of information 
and to extract correlations, ad-hoc algorithms and ap- 
propriate computer softwares have to be developed. Fur- 
thermore, once the space of the search is narrowed, a 
detailed study can be employed on the obtained reduced 
set of feasible candidate systems. 

Several examples have appeared in literature in re- 
cent years, for instance the "data-mining of quantum 
calculations" method leading to the principle component 
analysis of the formation energies of many alloys in sev- 
eral configurations 0, || , the evolutionary approach for 
determining hamiltonian Q, the "Pareto-optimal" al- 
loys and catalysts d, [^1 the prediction of the lithium- 
boron superconductor [lOj, the "high-throughput Kohn- 
anomalies" search in ternary lithium-borides [ll|, EH; 
and the "multi-optimization" techniques used in studying 
high-temperature reactions in multicomponent hydrides 

This manuscript focuses on the high-throughput for- 
malism for the calculations of solubility in binary alloys 
(solvus lines). The knowledge of solubility is crucial for 
designing new alloys with particular physical, chemical, 



and mechanical properties. For example, if we have to en- 
hance an alloy property by adding extra specie as solute, 
it is necessary to know the equilibrium solubility to un- 
derstand if it is possible to dissolve the candidate specie, 
and, if possible, to avoid supersaturation-precipitation 
and subsequent modification of the target property (ag- 
ing effect). In superconducting materials research, the 
problem emerges frequently: often expensive and diffi- 
cult experiments are undertaken to enhance the critical 
temperature [H, [lT} ■ Even in catalysis research, recent 
experiments and modeling have shown that the solubil- 
ity of carbon is responsible of thermodynamic instabili- 
ties hindering the catalytic activity of very small Fe and 
Fe:Mo clusters EE!]. The calculation of solubility of Zr 
in Al has already been addressed with success within the 
regular solution model (2fj ] fit to ah initio calculations 
[2lj . Solubility can also be extracted from the knowl- 
edge of the phase diagram which, in the case of lattice- 
conserving allo ys, can be generated within the the cluster 
expansion [22l |23{ and Monte Carlo approaches [24l [25|. 
However, a straightforward formalism leading to the esti- 
mation of equilibrium solubility for general alloys is still 
lacking. 

In the present paper we devised a statistical- 
thermodynamic approach for the calculation of atomic 
solubility in alloys. The advantage of our approach con- 
sists in taking into account all available ground states 
rather than just the pure species configurations. To 
test the method, we present calculations for a number 
of binary titanium systems. The paper is organized as 
follows: In Sec. II we rewrite the equations govern- 
ing solubility in the case of vacancies and substitutional 
impurities in binary alloys. Section III is devoted to 
the discussions of capabilities and limits of our formal- 
ism. Examples of phase diagrams and solubilities are 
addressed in Sec. IV for the following test Ti-A systems 
(A=Ag,Au,Cd,Co,Cr,Ir,W,Zn). Conclusions are given in 
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Sec. V. 



II. SOLUBILITY FORMALISM 
Enthalpy 

Let us consider a disordered dilute solution of A-atoms 
and vacancies (V) in a pure B-solid with a given Bravais 
crystal lattice (labeled with the "dis"). Without taking 
into account the A-A, A-V, and V-V interactions and 
assuming that A and V concentrations are small, the ap- 
proximate enthalpy of the considered alloy can be written 
as: 

H dis = E dis +pV dis = H *t N + h ^ Na + H ^ Nv: (1 ) 

where , H k ^ and i/v B are the enthalpy of the pure B 
solid per unit cell, the change in enthalpy of the solid 
upon substitution of one B with an A-atom, and the 
change in enthalpy upon removal of one B atom, respec- 
tively. In addition, N and N a (a=A,B,V) are the total 
numbers of crystal lattice sites and of atoms of a-type: 



N = N k + N B + Nv, = £g 



(2) 



B' 



and u° (a=A,B) are energy and the volume per atom 
of the a-pure solid, p is the pressure, v k and v2 rcprc- 

D B 

sent the change of volume of the B-purc solid upon intro- 
duction of one A-atom or one vacancy. The framework 
introduced by Eq. ([I]) is similar to Wagner-Schottky 
model of a system of non-interacting particles. [26j. The 
quantities with superscript "0" are considered to be tem- 
perature, pressure, and composition independent being 
calculated at zero temperature and pressure. Hence, the 
accuracy of the following results will be better in the case 
of limited pressures or for systems with very high bulk 
modulus, where the elastic energy fraction of is neg- 
ligible. 

The quantities in Eq. ([2]) can be easily approximated 
as differences of first-principles energies and volumes be- 
tween large supercells (sc) with or without defects (sub- 
stitutional A-atom or vacancy): 



E sc [B nsc ^ 1 A] - E sc [B nsc ] 
V sc [Bn sc -iA] - Vsc[-B nsc ], 



£sc[-Bn sc -l] - E SC [B nsc ], 

Vsc[B nsc -i] - V S c[-B„ sc ]. 



(3) 



As the size of the supercell grows, the approximate 
quantities in Eqs. Q approach their exact values. In 
literature and in this paper, E®, v® and H aR 



(a— A,V) are called the "raw" (composition unpreserv- 
ing) a-defect formation energy, volume and enthalpy, re- 
spectively. 

The enthalpy per atom is obtained from Eq. {!]) as: 



H dls /{N k + N B ) 



(4) 



+ H kji x k + (H Vb + H^)x v , 

where x a (a=A,B,V) are the atomic concentrations 

x a =N a /(N k + N B ). (5) 

Then, we follow the convention of using the formation 
enthalpy AH^ S calculated with respect to the pure A- 
and B-solids |28j : 

AH^ S = H£ s - x k Hf - (1 - x k )Hr , (6) 

where Hf* = Ef* +pv k and has been defined in Eq. 
(|2|). Combining equations (j4j and ([6]), we get 

AH^ S = H k x k + H v xv, (7) 

where the quantities H k and Hy, defined as 



B J 



Hv — H\r 



(8) 



are usually called A-defect and V-defect "true" (compo- 
sition preserving) [13 formation enthalpies, which can be 
obtained from [231 



dAH. 



dis 



dx a 



(a = A,V). 



(9) 



x a — >0 



Equilibrium Gibbs Free Energy 

The formation Gibbs free energy, AGat S ; is defined as: 

AGH S =AH£ S -TAS£ S , (10) 

where AH^ S is described by Eq. ([7|) and the formation 
entropy AS£t s = S^ s can be obtained within the mean- 
field approximation as 



AS'. 



dis 



k B N 
N k + N B 



c a In ( 



(11) 



a=A,B,V 



where T, k B , c a (a=A,B,V) are the temperature, the 
Boltzmann constant, and the site concentrations of 
atoms: 

c a = NJN. (12) 

By changing variables from site concentrations c a to 
atomic concentrations x a : 

x a = N a I (N k + N B ) , /(1 + xv), (13) 
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we rewrite Eq. (|10|) as 



AG 
AG' 



dis, A 
at 

dis,V 
at 



H^Xk + k B T[x k lnx A + (1 - xa) hi(l - x A )], 



iJyxv + kaT [xt/ In xy — (1 + Xv) ln(l + xv)]. 

(14) 

In alloy with fixed atomic composition xa, the equi- 
librium concentration of vacancies x^ is determined by 
minimizing the formation Gibbs free energy: 



dAG 



dis 



dxy 



dAG 



dis.V 



'A 



dxy 



0. 



The manipulation of Eqs. (fT¥|) - (fT5|) leads to: 



„ e< l 



„ ec i| 



exp 



1 — exp 



(15) 



(16) 



AG 



i dis g S i 




/ X 1 

/ dis+gs V 




4 H 




J i dis gs i 


T=T o 


v \ 


K 


1 '*"* • \ 1 
i ■ -« . \ i 




(b) 



B 



x 501 x 

A *A 



AH, 



agS"' t (^) 



~ -£; B Texp ( 

fe B T<B v V B 



To conclude, the "true" vacancies formation enthalpy 
flv determines the equilibrium concentration of vacan- 
cies with an Arrhenius-type equation (see also Ref. l28h . 
In the next section, we show that resolves the solu- 
bility in the case of a phase-separating alloy having no 
intermediate ground states. 



Solubility 

At a given temperature, the solubility of A-atoms in a 
B-solid phase, xf 01 , is defined as the maximum homogc- 
nously achievable concentration of A without the forma- 
tion of a new phase (Fig. QJa)). The accurate calculation 
of x| o1 requires considering of the nearby ground state 
(labeled as "gs") with respect to the increase of xa, see 
the "red line" in Fig. [TJb). It is implemented by min- 
imizing the formation Gibbs free energy AG^ x (x) of a 
mixture of (1) a disordered dilute solution of A-atoms 
and vacancies in a B-rich solid phase at composition xa b 
(the "dis" -phase of the previous section), and (2) an on- 
or off-stoichiometric ground state "gs" -phase at compo- 
sition XAg S • The lever rule gives the fractions of the two 
phases: 



AG^(x) = 



^AGK(* Ags ), 



(17) 



FIG. 1: (Color online) (a) Temperature-concentration and (b) 
formation Gibbs free energy-concentration (at a given temper- 
ature To) graphs illustrating our solubility concepts. 



and the minimization is performed with respect to xa b 
and XAgs (x is the overall composition of A in the two- 
phase mixture, xa b < x < XA gs ): 



<9AG^ X 
<9xa k 



0, 



<9XA e 



(18) 



Combining Eqs. (fTT)) and (fT5)) leads to the usual 
common-tangent rule: 



AG 



B 

dis 
at 



dx t 



(*A c )-AGg>Ags_) 



(19) 



Ag A gs 

Substituting Eqs. (Tl"4"|) into Eqs. p^|) . we obtain 
xff = [exp (JW/W + l]" 1 , 



(20) 



which approximates as an Arrhenius-type relation at low 
temperature: 



xif 



exp(-iJ so i/fc B r) 



(21) 
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The quantity H so i is defined as 



H so i = i?A — 



„,sol „,sol 
x k B x A gs 



(22) 



The non-linear problem described by Eqs. (|20II21|) and 
(|2"2"j) can be linearized in the low-solubility limit (labeled 
as "Is"): 



,,sol 



0, AGH*(xl?)~AGil s '\x«% 



r sol 



AGf t s (4°i) ~ Atff t s , 



(23) 



where A_ff^ is the formation enthalpy of the ground 
state "gs". In the low-solubility limit, H so i becomes: 



ls.nv 



^sol — -^so] 



where H^oi^ is the non- vacancy contribution (labeled as 



AG 



dis.V/ eq 



at 



(24) 



H 



ls,nv 
sol 



H A - AH^/i 



igs- 



(25) 



From the equilibrium vacancy concentration, Eq. (|16[). 
the exponential part of Eqs. (|20j|2ip becomes: 



cxp (-HlSJkBT) = exp(-H^/k B T) x 



1 — exp 



3l 



(26) 



where the two contributions, non- vacancy and vacancy, 
are factorized. The last expression indicates that the 
presence of vacancies effectively increases the solubility 
by decreasing the number of host B-atoms in the solution. 



III. INTERPRETATION OF H 



ls,nv 
sol 



For low-solubility calculations of non-interacting de- 
fects, the framework can be implemented through first- 
principles calculation of Hl^i (-^ < 4- p4|) )- It requires the 
knowledge of the enthalpy Ai/f{° and composition xl 
of the ground-state "gs" , as well as the "true" defect for- 
mation enthalpies H& and Hy (Eq. ([5])). 

To capture the physical meaning of H^^ v , let us con- 
sider an arbitrary dilute solution "sc" at composition xf c 
(see Fig. [2]). The label "sc" indicates that the solution is 
generated as a supercell of the B-solid upon insertion of 
defects, randomly distributed but not too close to avoid 
interactions. The first part of Eq. (|2"5|) can be rewritten 
according to Eq. ([9]) as: 



dxf c 



(27) 




-AH 



AH SC -A/-/ ( (B+9S) 
at at 



gs 



H 



H 



(B+gs) 



at 



B 



FIG. 2: (Color online) The enthalpy-concentration graph 
demonstrating our concept of low-solubility formation en- 
thalpy. A-pure, B-pure and "gs" are the ground states form- 
ing a convex hull, "sc" is a compound (supercell) with some 
intermediate enthalpy and composition. Compounds (B+A) 
and (B+gs) correspond to appropriate phase mixtures of the 
same general composition a;| c . 



The second part of Eq. (|25[) becomes: 



Afff t s 



ftgs 



AH. 



( B +g S )^SC^ 



dAH. 



(B+gs; 



dx\ c 



(28) 



where Ai7^ +gs) (xf 



is the formation enthalpy of the 
mixture of the pure B-solid with the ground state "gs" 
with overall composition xf c (point (B+"gs") in Fig. [2]). 
Thus, we obtain: 



H. 



ls,nv 
sol 



d 



AHg(xi c )-AHg +ss \xl c ) 



(29) 



The comparison of Eq. (g7| with Eq. (gSJ) leads to the 
conclusion that both H& and H^^" are derivatives of su- 
percell formation energies with respect to A-composition. 
For Ha, the supercell formation enthalpy is determined 
with respect to pure A and B phases (the distance be- 
tween points "sc" and (B+A) in Fig. [2]). For f/ 'j™, the 
supercell formation enthalpy is determined with respect 
to B-pure and the ground state "gs" (the distance be- 
tween points "sc" and (B+gs) in Fig. [2]). To conclude, 
and H^^™ are characterized by the angles a and 
(3 between the B-"sc"/B-A and B- "sc" /B- "gs" lines, re- 
spectively. In analogy with the definition in Eq. ([S]), 



H 



ls.nv 
sol 



can be considered as the "low-solubility formation 



enthalpy" . 

The quantities Hp, and H^^" are identical only for 
phase-separating alloys having no intermediate ground 
states (i.e. "gs"= A). In this case the low-solubility 
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can be formally determined by minimizing AG^t S ' A (Eq. 
(fl4|) ) with respect to x& (e.g. Ref. HH). However, in the 
general case, the existence of ordered ground-states must 
be verified so the appropriate formalism is used. Gener- 
ally, in ordering alloys and H^^ v differ, and might 
even have different signs. 

Note that if the "sc" point is below (B+gs) in Fig. [2] 
(#soi nV < °) the solubility expressions (|20H21| are not 
valid, and there must exist an undetected ground state 
(it might be "sc" itself) with concentration lower than 
%1 SS . In this case, such undetected ground state should 
be used for the calculation of solubility rather than initial 
"gs" & 

The expression for non-binary low-solubility within the 
regular solution model derived in Refs. [2(| El coincides 
with our derivation in the case of binary alloys without 
vacancies and high-temperature contributions. This is 
because the regular solution model corresponds to our 
model for the free energy in case of dilute solution. 



IV. GROUND STATES OF SELECTED 
TITANIUM ALLOYS 

As an example of our formalism, we calculate the sol- 
ubility of a set of metals in titanium. First, we ex- 
plore the possible ground states of the A-Ti systems 
(A=Ag,Au,Cd,Co,Cr,Ir,W,Zn) and then we apply the 
construction described in the previous section. 

The low temperature stability of A-Ti is performed by 
using our hi gh-t hroughput quantum calculations frame- 
work 0, [H, lid [HJ , based on first-principles energies 
obtained with the VASP software [3l|. We use pro- 
jector augmented waves (PAW) pseudopotentials HH 
and exchange-correlation functionals as parameterized by 
Perdew-Burke-Ernzerhof [33| for the generalized gradient 
approximation (GGA). Simulations are carried out with 
spin polarization except Ti-Co, at zero temperature and 
pressure [34| . and without zero-point motion. All struc- 
tures are fully relaxed (shape and volume of the cell and 
internal positions of the atoms). The effect of lattice 
vibrations is omitted. Numerical convergence to within 
about 1 me V/ atom is ensured by enforcing a high energy 
cut-off (357 eV) and dense 6,000 k-point meshes. 

The number of crystal structures considered for the cal- 
culations of each A-Ti system is 194. In addition to the 
176 described in Ref. 0, we included the following proto- 
types A5, A6, A7, All, Ca 7 Ge, NbNi 8 (Pt 8 Ti), V 4 Zn 5 , 
C36 and the whole complete set of hep-superstructures 
with up to 4 atoms/cell. The whole process is performed 
in an automatic fashion through the software AFLOW 
which generates the prototypes, optimize the parameters, 
perform the calculations, correct possible errors, and cal- 
culate the phase diagrams 0, [35[ . 

Our results of ground state calculation are presented 
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Au Atomic Percent Titanium Ti 

FIG. 4: (Color online) AuTi (gold-titanium) ground state 
convex hull. 



in Figs. [BTfTfJl and Table fl] The correspondence ab initio 
versus experimental results is very good and typical for 
this type of calculations [j| . 

We propose degenerate results for: 1) Co2Ti, experi- 
mentally reported as C15, but with ab initio formation 
energies of -311.4 meV/at. and -304.0 meV/at. for C14 
and C15, respectively; 2) CoTi2, experimentally reported 
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FIG. 5: (Color online) CdTi (cadmium-titanium) ground 

state convex hull. 
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FIG. 6: (Color online) CoTi (cobalt-titanium) ground state 

convex hull. 



0.02 




FIG. 7: (Color online) CrTi (cromium-titanium) ground state 
convex hull. 
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Ti Atomic Percent Tungstem W 

FIG. 9: (Color online) TiW (titanium-tungsten) ground state 

convex hull. 
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FIG. 10: (Color online) TiZn (titanium-zinc) ground state 
convex hull. 



as NiTi2, but with ab initio formation energies of -291.2 
meV/at., -287.3 meV/at., and -285.7 meV/at. for C37, 
CuZr2, and NiTi2, respectively; 3) TiZn, experimentally 
reported as B2, but with ab initio formation energies of 
-195.4 meV/at. and 193.4 meV/at. for Ll and B2, re- 
spectively. 

We propose the novel results for: 1) IryTi, experi- 
mentally reported as a two phases region above 500° C, 




20 40 60 80 100 

lr Atomic Percent Titanium Ti 

FIG. 8: (Color online) IrTi (iridium-titanium) ground state 
convex hull. 



but with a possible ab inito low temperature ground 
state CayGe with formation energy of -373.8 meV/at. 
2) Ir2Ti, experimentally reported as a two phases re- 
gion above 500°C, but with a possible ab initio low tem- 
perature ground state C11&, with formation energy of 
-716.0 meV/at. 3) T1W2, experimentally reported as a 
two phases region above 500°C, but with a possible ab 
initio low temperature ground state BCCj^ 1 ^, with for- 
mation energy of -82.7 meV/at. (see notation for the 
prototype in Ref. [f|). 4) T1W4, experimentally reported 
as a two phases region above 500° C, but with a possible 
ab initio low temperature ground state Dl a , with forma- 
tion energy of -83.8 mcV/at. 5) Ti3Zn, experimentally 
not explored, but with a possible ab initio low tempera- 
ture ground state A15, with formation energy of -120.0 
meV/at. 

In particular, the results indicate that in the Ir-Ti sys- 
tem the low temperature Ir-rich part of the known phase 
diagram is not complete, and that the Ti-W alloy has 
an ordering tendency at low temperature, in contrast to 
common belief (36l.[37j. In addition, in Ti-Zn there must 
exist a Ti-rich compound with Ti composition higher 
than the reported T^Zn, [3H 
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TABLE I: Low temperature phases comparison chart for Ag- 
Ti, Au-Ti, Cd-Ti, Co-Ti, Cr-Ti, Ir-Ti, Ti-V, and Ti-Zn. Ex- 
perimental data correspond to the lowest available tempera- 
ture 0| . The "*" indicates the energy values used in solubility 
calculations. 



Experimental 


Ab initio 




Space 


Pearson 


(Refs. 36, 37) 


result 


(meV/at.) 


group [38] 






Ag 


I-Ti 






AgTi-Bll 


Bll 


-67.6 


P4/nmm 


tP4 


AgTia-ClH 


Cllfc 


-63.3* 


I4/mmm 


tI6 


Au-Ti 


Au 4 Ti-Dl a 


Dl n 


-283.2 


I4/m 


tno 


Au 2 Ti-Cll b 


Clli 


-430.4 


I4/mmm 


tI6 


two-phase region 


Au 4 Ti 3 - 


-430.6 


I4/mmm 


tI14 


above 500° C 


Cu 4 Ti 3 /tic 








AuTi-Bll 


Bll 


-429.8 


P4/nmm 


tP4 


AuTi 3 -A15 


A15 


-356.1* 


Pm3n 


cP8 


Cd-Ti 


CdTi-Bll 


Bll 


-62.0 


P4/nmm 


tP4 


CdTia-Clli, 


Cllj 


-73.7* 


I4/mmm 


tI6 


Cr-Ti 


Cr 2 Ti-C15 


C15 


-120.3* 


Fd3m 


cF24 


Co-Ti 


C03T1-LI2 


Ll 2 


-254.8 


Pm3m 


cP4 


Co 2 Ti-C15 


C14 


-311.4 


P6 3 /mm 


hP12 




C15 


-304.0 


Fd3m 


cF24 


CoTi-B2 


B2 


-386.4 


Pm3m 


cP2 


CoTi 2 -NiTi 2 


C37 


-291.2* 


Pnma 


oP12 




CuZr 2 


-287.3 


I4/mmm 


tI6 




NiTi 2 


-285.7 


Fd3m 


cF96 


Ir-Ti 


two-phase region 


Ir 7 Ti- 


-373.8 


Fm3m 


cF32 


above 500° C 


Ca 7 Ge 








Ir 3 Ti-Ll 2 


Ll 2 


-720.3 


Pm3m 


cP4 


IrTi-<5 


Llo 


-851.8 


P4/mmm 


tP4 


two-phase region 


Ir 2 Ti- 


-716.0 


I4/mmm 


tI6 


above 500° C 


Clli 








IrTi 3 -A15 


A15/tie 


-570.3* 


Pm3n 


cP8 


Ti-W 


two-phase region 


TiW 2 - 


-82.7* 


P3ml 


hP3 


above 500° C 


Rrr ,[211] 








two-phase region 


TiW 4 - 


-83.8 


I4/m 


tno 


above 500° C 


Dl n 








Ti-Zn 


non explored 


Ti 3 Zn-A15 


-120.0* 


Pm3n 


cP8 


Ti 2 Zn-CuZr 2 


CuZr 2 


-158.0 


I4/mmm 


tI6 


TiZn-B2 


Llo 


-195.4 


P4/mmm 


tP4 




B2 


-193.4 


Pm3m 


cP2 


TiZn 2 -C14 


C14/tie 


-197.2 


P6 3 /mm 


hP12 


TiZn 3 -Ll 2 


H a 


-198.0 


Pm3m 


cP4 



V. RESULTS AND DISCUSSION OF 
SOLUBILITY IN TI ALLOYS 

The "raw" formation enthalpies of Ti alloys with one 
substitutional atom or a vacancy were obtained through 
Eqs. ([3]) considering 3x3x3 supercells of the hep Ti. Su- 



percell dimensions were chosen to limit the defect-defect 
interactions, being their distance at least three times 
larger than the nearest neighbor Ti-Ti bond. The re- 
sults are presented in Table [UJ Solubilities temperature 
dependencies are presented in Fig. [Til while Fig. [12] 
illustrates a comparison of experimental and theoretical 
data at T=700°C. 

TABLE II: "Raw" formation enthalpies , ^v Ti (from 

Eqs. d2][3j) at zero pressure), "true" formation enthalpies 
Hy (from Eq. at zero pressure) and low-solubility (non- 
vacancy) formation enthalpy H ^ v (from Eq. ([25} at zero 
pressure) of substitutional A-defects (A=Ag,Au,Cd,Co,Cr) 
and vacancies (V). The elements are ordered from low to high 
"low-solubility formation enthalpy". All quantities are in eV 
units. 



A 




H A 


TT-Ls.nv 


Zn 


6.405 


-0.262 


0.218 


Cd 


7.267 


0.24 


0.461 


Ag 


5.412 


0.304 


0.494 


Au 


3.882 


-0.781 


0.643 


W 


-4.303 


0.714 


0.838 


Ir 


-2.059 


-1.14 


1.141 


Co 


1.138 


0.316 


1.19 


Cr 


-0.674 


1.024 


1.205 


H Vt . = 10.002, 


#V = 


2.068 



Table[TTJshows negative "true" formation enthalpies if a 
for A=Zn, Au, and Ir. This indicates that calculation of 
Zn, Au, and Ir solubilities in Ti is not faceable without 
considering nearby intermctallic ground states. 

Figures [TTTfTSl show that the highest theoretical solu- 
bilities in Ti occurs for Zn, Cd, Ag, and Au as conse- 
quence of their low solubility formation enthalpies (see 
Table [TlJ . This high solubility has also been observed 
experimentally for Cd, Ag, and Au, whereas, to our best 
knowledge, solubility of Zn does not seem to have been 
studied. High solubility of late transition metals Zn, Cd, 
Ag, and Au in the early transition metal titanium is due 
to the substantial localized stability provided by the fill- 
ing tendency of the d-band of Ti (Ref. l39h . 

Because of the very high formation enthalpy of vacan- 
cies in Ti (reported in Table Hi*)) . the vacancy equilibrium 
concentration was found to be very low at all consid- 
ered temperatures (e.g. < 1(T 6 at T < 1300°C). 
Thus, the effect of vacancies on the solubilities of Ti al- 
loys should be negligible. 

Although, theoretical and experimental results follow 
similar trends, theoretical solubilities are considerably 
smaller for most of the considered alloys. A similar dis- 
crepancy was also observed for Al-Zr in Ref. [2l|. The dis- 
crepancy could be due to shortcomings of theory and/or 
experiment. The main approximations of our model con- 
sist of (a) neglecting the defect interaction, (b) neglect- 
ing the spatial defect correlation and (c) assuming low 
concentration of defects. However, as such approxima- 
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T(°C) 

FIG. 11: (Color online) The solubilities of eight considered 
transition metals in titanium as functions of temperature (cal- 
culated through Eqs. (|2QI24p h Panel (b) is a magnified ver- 
sion of panel (a), to visualize the values for Ir, Co, and Cr. 



tions are somehow related, subsequent solubility calcula- 
tions suggesting low values validate the assumptions (the 
mean-field approximation neglects the interatomic posi- 
tional correlations but it should work well when the devi- 
ation from complete stoichiometric (pure Ti-solid in our 
case) is small - and so is the solubility in our case - see Sec. 
19 in Ref. H^.) For all considered alloys except Ti-Zn, we 
found that at intermediate temperatures the calculated 
solubilities arc small enough for our approximations to be 
valid, although smaller than experiment, as mentioned 
before. In RefiH, the authors added the defect interac- 
tions through the Cluster Expansion method but with- 
out increasing the solubility considerably Actually, our 
formalism includes the solute-solvent ordering tendency 
by considering the real intcrmctallic ground state other 
then the pure Ti-solid. Thus, we conclude that our ap- 
proximations (or those of the cluster expansion method) 
are not responsible for the theory-experiment solubility 
discrepancy. 

In case of Ti-Zn, the formation enthalpy is very low 
(see Table [TTJ) . Correspondingly, the theoretically solu- 



5 
4 

o 

>3 



r=700°c 

O— Experiment 
•—Theory 




Zn 



Cd Ag 



Au 



W 



Ir 



Co Cr 



FIG. 12: (Color online) Comparison of experimental [39. 137] 
and theoretical solubilities of eight considered transition met- 
als in titanium at T=700°C. From left to right, the elements 
are ordered from low to high low-solubility formation enthalpy 
(and correspondingly theoretical solubility) The error bars 
characterize the scattering of data measured in different ex- 
periments. We could not find the experimental data for solu- 
bilities in Ti-Zn. 



bility at intermediate temperatures is very high, violat- 
ing the assumptions of the model. Thus, for Zn in Ti, a 
more precise solute interaction and correlation parame- 
tcrizations is required, which can be obtained by using, 
for example, Cluster Expansion and Monte Carlo simu- 
lations. Hence, the presented formalism should be con- 
sidered as an appropriate starting point to determine if 
more computationally expensive formalisms are needed. 

The other approximation of our model is the assumed 
independence of our model energy and volume parame- 
ters on temperature. The dependence can be caused by 
non configurational degrees of freedom, like vibrational 
or anharmonicity (magnetic ordering is not actual for 
considered alloys). However, theory-experiment solubil- 
ity discrepancy is observed even at low enough temper- 
atures (e.g. T < 9 D (Ti) =374-385 K [H|), wh ere the 
vibrational contribution to the free energy is not impor- 
tant. In fact, vibrational entropy is substantially smaller 
than the configurational contribution [42j , so its inclusion 
can not modify solubility results much. In fact, even in 
Ref. I2H the authors added phonon contribution without 
increasing the solubility considerably. 

On the other side, the experimental equilibrium solu- 
bility tends usually to be overestimated. In fact, the for- 
mation of metastablc and/or unstable states which are 
subsequently frozen at low temperatures, can make solu- 
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bility measurements very challenging. In such scenarios, 
the measured solubility may correspond to spinodal con- 
centration rather then actual binodal concentration or 
simply characterize the frozen out of equilibrium solubil- 
ity remaining from the initial specimen preparation at 
higher temperature. Besides, the segregation of defects 
into grain-boundaries, especially in multicrystalline sam- 
ples prepared through non optimal cooling dramatically 
affect the amount of frozen defects and solutes. 



VI. CONCLUSIONS 

Based on the statistical-thermodynamic theory of di- 
lute lattice gas, we developed an approach to calculate 
the atomic bulk solubility in alloys. The advantage of 
our approach consists in considering all the appropriate 
ground states rather than the pure species. It was shown 
that the low-solubility follows simple Arrhcnius-type 
temperature dependence determined by a "low-solubility 
formation enthalpy". This quantity is defined as the com- 
position derivative of the compound formation enthalpy 
with respect to nearby ground states. "Low-solubility 
formation enthalpy" coincides with the usual defect for- 
mation enthalpy only in the case of a phase-separating 
alloy having no intermediate ground states and vacancies. 
The key quantities of our model can be directly obtained 



by first-principles calculations of by fitting experimen- 
tal temperature solubility dependence. Generalization of 
our model to intermediate phases and/or to multicompo- 
nent, multisublattice, interstitial-substitutional alloys is 
straightforward. 

As examples, we applied the framework for a set ot 
eight Ti alloys A-Ti (A=Ag,Au,Cd,Co,Cr,Ir,W,Zn). We 
have found that the highest solubility for Zn, Cd, Ag, 
and Au is in qualitative agreement with available ex- 
perimental data and band structure expectations. The 
quantitative differences between the theory and experi- 
ment observed in the present and other similar studies 
are discussed. 

In conclusion, our formalism is correct in the limit of 
low-solubility and should be considered as an appropriate 
starting point for determining if more computationally 
expensive formalisms are otherwise needed. 
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